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PREFACE 


In 1955, the team which has become the Marshall Space Flight Center 
(MSFC) began to organize a research program within its various laboratories 
and offices. The purpose of the program was two-fold: first, to support existing 
development projects by research studies and second, to prepare future develop- 
ment projects by advancing the state of the art of rockets and space flight. 
Funding for this program came from the Army, Air Force, and Advanced 
Research Projects Agency. The effort during the first year was modest and 
involved relatively few tasks. The communication of results was, therefore, 
comparatively easy. 

Today, more than ten years later, the two-fold purpose of MSFC’s 
research program remains unchanged, although funding now comes from NASA 
Program Offices. 'The present yearly effort represents major amounts of money 
and hundreds of tasks. The greater portion of the money goes to industry and 
universities for research contracts. However, a substantial research effort is 
conducted in house at the Marshall Center by all of the laboratories. The com- 
munication of the results from this impressive research program has become a 
serious problem by virtue of its very voluminous technical and scientific content. 

The Research Projects Laboratory, which is the group responsible for 
management of the consolidated research program for the Center, initiated a 
plan to give better visibility to the achievements of research at Marshall in a 
form that would be more readily usable by specialists, by systems engineers, 
and by NASA Program Offices for management purposes. 

This plan has taken the form of frequent Research Achievements Reviews, 
with each review covering one or two fields of research. These verbal reviews 
are documented in the Research Achievements Review Series. 


Ernst Stuhiinger 

Director, Research Projects Laboratory 


These papers presented September 16. 1965 
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INTRODUCTION TO RESEARCH ACHIEVEMENTS 
REVIEW ON AERODYNAMICS RESEARCH 


By 

Dr. E. D. Geissler* 


The subjects in this aerodynamics review are 
varied and cover such topics as prediction methods 
of drag and lift characteristics, jet penetration in 
water, base -flow phono mena, optimum nose shape 
design and several topics related to aerodynamic 
measurements. The paper by Dr. F. Krause dis- 
cusses a very promising new optical cross-correlation 
method. If the expected results of this method ma- 
terialize, the state-of-the-art of some aerodynamic 
measurement techniques will be advanced markedly. 

Dr. Platzer's contribution concerns some major 
research activities in unsteady aerodynamics, a 
field which represents still several ticklish problems 
for the Saturn program (in particular, aerodynamic 
noise-and-ground-wind oscillation) . Mr. Homer 
Wilson reports on aerothermodynamic research 
studies, primarily base heating, which is not yet in 
a very satisfactory shape in spite of substantial ex- 
perimental efforts. This is particularly true when 
it is necessary to predict new configurations. 

A final paper by Mr. Jim Ballance gives a short 
discussion of studies on rarefied-gas dynamics - a 
topic which, while not yet so prominent in producing 
engineering headaches for the designers of space 
vehicles, will undoubtedly be of growing concern with 
extended stay time in very high altitudes and which 
especially poses many unsolved questions concerning 
instrumentation. 

The activities reported are all carried out with 
participation and supervision of personnel of the 
Aerodynamics Division, even though a major share 

* Director, Aero-Astrodynamics Laboratory. 


is done by contractors in Government and industry. 

It is hoped that these presentations will succeed in 
impressing the reader with the scope and importance 
of aerodynamics research in space-vehicle design 
and will also illustrate the particular flavor of applied 
aerodynamic research which is in close connection 
with the development of space vehicles. Problems 
are stimulated by the need for practical answers and 
sometimes by results from flight tests rather than 
selected according to their academic appeal. The 
justification of aerodynamic research facilities at the 
vehicle development agency, besides the larger aero- 
dynamic research center facilities, stems from 
similar reasons and appears amply demonstrated by 
efficiency in tackling immediate problems. 

Aerodynamics of space vehicles are relatively 
ungl amorous compared with, for example, develop- 
ment of power plants. With increasing size of ve- 
hicles, the immediate importance of aerodynamic 
forces versus inertial effects decreases. Whereas 
the performance of an airplane is the result of a 
direct balance between aerodynamic drag and lift 
characteristics versus available thrust, the perform- 
ance for a big space vehicle is mainly defined by 
propulsive inertial balance. However, it must be 
kept in mind that the efficiency of the structural de- 
sign causes increasing concern with increasing rocket 
size and this structural design is strongly influenced 
by aerodynamic loads (i. e. , bending loads, necessary 
control forces, aerodynamic noise, base heating). 

The efficient design frequently favors shapes which 
are not conducive to undisturbed flow characteristics 
(protuberances and separation including geometry) 
and requires prediction of flow not encountered in 
airplane work. 



GENERAL REVIEW OF AERODYNAMICS RESEARCH 


By 

Werner Dahm 


I. INTRODUCTION 


The scope of aerodynamics research at MSFC is 
determined by the fact that, in accordance with the 
Center’s mission, it has to serve space -vehicle de- 
velopment, and especially launch-vehicle development. 
Under this guideline, this work includes: (a) de- 

velopment or improvement of analytical methods, and 
their adaptation to our needs; (b) experimental in- 
vestigations; (c) development and improvement of 
experimental methods and facilities; (d) development 
of instrumentation and probing tools; and (e) general 
design studies. 

As a further guideline, tasks which do not 
promise useful design applications within a period of 
about 5 years are considered beyond MSFC f s scope. 
The majority of the tasks yield results which are 
applied within a much shorter period. 

The following are some typical examples of our 
research work; some other areas will be covered in 
separate papers. 

(a) Subsonic flow. 

(b) Jet impact in water. 

(c) Base-flow investigation. 

(d) Effect of wall-to-total temperature ratio on 
hypersonic flow detachment. 

(e) Variable porosity walls. 

(f) Dynamic balance for Saturn forebody. 

II. SUBSONIC FLOW 

The aerodynamic characteristics of the Jupiter 
nose cone which are shown in Figure 1 are examples 
of the analytical work. This example is concerned 
with pressure and lift distribution on bodies of re- 
volution in subsonic flow. The bodies which are 


dealt with always have a blunt base, a case which 
available theories could not properly handle. There- 
fore, a method was developed which consists of dis- 
tributing singularities on the surface of the body, and 
using the special boundary condition that the rim-base 
pressure is equal to the base pressure, and that the 
local lift is zero. The base pressure cannot be given 
by inviscid theory, and must be known a priori; the 
zero-lift condition is an equivalent of the kutta- 
condition of wing theory, and reflects the fact that 
virtually no pressure difference can be maintained 
across the ’’dead air” of the wake. 

Figure 1 compares analytical and experimental 
results on a blunt-nosed cone. A blunt-nosed cone 
was selected because a complete set of experimental 
data were available, and because the base-effects are 
fairly pronounced. The zero -lift pressure distribution 
in the lower left corner coincides with the experimental 
results of Figure 1. This is also true for the lift dis- 
tribution in the upper left corner of Figure 1. At this 
particular Mach number, lift and pressure over the 
rear 20 percent of the body are noticeably influenced 
by the base. As the Mach number decreases, the base 
influence extends further forward, and lift is lost on 
the rear of the body. The corresponding effects on 
the total lift and center of pressure can be seen on the 
2 graphs on the right of Figure 1. Since the method 
correctly describes the flow field near the base in- 
cluding the effects of the base pressure, the effects 
of the rocket jets on the vehicle stability can be es- 
timated once the base pressure with jets on is known. 

Figure 2 shows another example of the application 
of this method in that it depicts the drag of cones with 
tip half angles from 15° to 90® as function of the base 
pressure. The experimental points were obtained 
from a water tunnel, where the base pressure is 
fixed by the cavitation conditions; therefore, the base 
pressure coefficient can be varied. Agreement be- 
tween analyses and experimental results is very good. 
Closer inspection shows that the total drag varies less 
than the base drag, reflecting again the forward effect 
of the base pressure. Incidentally, the analysis yields 
a correct description of the major part of the cavi- 
tation bubble. This result is of interest to submarine 
applications; to us, it is only an interesting by-product. 
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III. JET IMPACT IN WATER 


Jet penetration into water represents another 
example of our analytical work (Fig. 3). There are 
a number of reasons why launch facilities should be 
placed in the offshore waters such as launch-site 
real estate requirements. In this case the rocket 
jets would be discharged into the water, and it is 
necessary to know how deeply they penetrate. Sc ale - 
model simulation is practically impossible (it would 
require testing in a low-pressure chamber, or on a 
rocket sled under controlled, high acceleration) . 
Therefore, the question has been attacked by com- 
posing a simple, approximate model of the process, 
and by checking the. validity of this model with the 
results of small-scale jet tests. The model which 
was used was the well-known variation of the impact 
pressure of turbulent free jets with the distance down- 
stream of the nozzle exit, and the fact that: (a) the 

penetration depth must correspond to the local impact 


M * 0 



| £■ Disk p - 90 


Cone p * 45 


Cone p • 15° 


FIGURE 2. DRAG COEFFICIENT OF CONES OF 
DIFFERENT APEX AS A FUNCTION OF THE BASE 
PRESSURE COEFFICIENT 
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IV. BASE FLOW INVESTIGATION 



FIGURE 3. PENETRATION DEPTH OF SUB AND 
SUPERSONIC JETS IN WATER 


pressure of the jet; and (b) the buoyancy of the 
cavity volume should equal the thrust of the engine. 
Figure 3 shows the correlation of a sizeable number 
of experiments with hot and cold subsonic and super- 
sonic jets, using the resulting relations. The abscissa 
represents the distance of the apex of the cavity from 
the nozzle exit. The ordinate contains the thrust F, 
the specific weight of the liquid y , the jet diameter d, 
the length of the in viscid jet core x^, the penetration 

depth n , and a Mach number function. The horizontal 

lines correspond to cases where the inviscid core of 
the jet hits the bottom of the cavity. The inclined 
line corresponds to cases where H + n^ > x^, i. e. , 

where the ’’inviscid core” of the jet no longer reaches 
the bottom of the cavity. The correlation of the data 
in Figure 3 shows inherently that the presence of the 
water has essentially no effect on the impact pressure 
profile of the jet. The meaning of these data may be 
illustrated by an application to the F-l engine. 


H 

n 

0 

r x engine. — 

0 

100'= 30.4 m 

160* = 50 m 
93’ = 28 m 


For closely clustered jets, the penetration depth 
is somewhat larger than for single jets and can be 
estimated from a model test of the impact pressure 
profile in free air. 


The various base-heating problems encountered 
in high speed reentry as well as from our engine jets 
have, for a long time, aroused our curiosity about 
what really goes on in the so-called ”dead-air zone” 
at the base. These processes have been studied on 
a two-dimensional, blunt based wing model at M = 3 
with a turbulent boundary layer. Figure 4 shows a 
shadowgraph picture of the base area of this model. 

The white rectangle on the left is the rear end of the 
model, 15 inches (37. 1 cm) high. The traces of the 
boundary layer, the Prandt- Meyer fan, the free shear 
layers, the tracking shocks, etc. can also be seen. 

Heat transfer and pressure distribution tests have 
been made on the base plate. The recirculation field, 
the free shear layers, and the base plate boundary 
layer have been probed with total pressure and tem- 
perature probes, wind vanes, an interferometer, and 
hot-wire anemometers. The lower part of Figure 4 
shows the density and distribution at zero heat flux, 
taken from an interferometer test. The flow in the 
recirculation zone was found to be extremely turbulent, 
/ 2 

withvu 1 = 0(u). The wind-vane tests of the mean 
velocity directions, shown by the arrows in the upper 
parts of Figure 4, deteriorated near the model center- 
line, indicating large outward flow angles rather than 
the required direction parallel to the centerline. The 
presumed cause is the high level of turbulence com- 
bined with strong shear flow effects. The main re- 
sistance to the heat transfer was found vested in a 
very thin base plate boundary layer, about 0. 06 inches 
(0. 152 cm) thick; the free shear layers played only a 
minor role. Outside of this base boundary layer, the 
total temperature is constant throughout the recircu- 
lation zone and the inner part of the shear layers. The 
base boundary layer is unlike any ordinary boundary 
layer. Its temperature profile seems to consist es- 
sentially of two straight lines, an outer one of small 
gradient dt/dx, and an inner one with steep gradient, 
which corresponds to the molecular heat transfer and 
reflects the existence of a laminar sublayer whose 
thickness is 0. 005 inch (0. 0127 cm) . The cause of 
this peculiar profile is believed to be the extreme 
turbulence level of the recirculation zone. Figure 5 
shows the effect of the base plate temperature on T^, 
the total temperature in the recirculation zone. T^ 

is seen to be very little affected by the base plate 

temperature T . dT/dT indicates the distribution 
1 w b w 

of the resistance to the heat flux between the free 
shear layers and the base boundary layer. 
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Recirculation Zone Temperature versus Base Plate Temperature 


FIGURE 5. TWODIMENSIONAL BASE FLOW 


V. EFFECT OF WALL-T0-T0TAL TEMPERATURE 
RATIO ON HYPERSONIC FLOW 
DETACHMENT 


To assess factors that have to be observed for 
proper wind-tunnel simulation, the wall-to-total tem- 
perature ratio was varied in a wind-tunnel test of a 
Saturn forebody model at a Mach number of 6. Figure 

6 shows schlieren pictures of the model at T /T 

wo 

- 0. 88 and 0. 22, demonstrating differences of flow 
detachment on the escape rocket-vehicle nose com- 
bination. Figure 7 shows the effect of T /T on the 

wo 

drag coefficient. At T /T near 1, the flow is fully 
wo 

detached and the drag is low. As the wall temperature 
decreases, the flow detachment recedes, increasing 
C^, until the flow is fully attached and has reached 
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Spike Shock' 


M = 6. 00 ; a * 0* M - 6. 00 ; a = 0* 

OO oo 

Schlierens of SATURN IB/Apollo/2nd Stage Flight Configuration 
s Weak 

Reattachment —j ron £ Reattachment I- 

Bow 

Spike Shock - 


Laminar 

Separation 


Bow 


Laminar 

Separation 


Engineering Interpretations of Schlierens 

Effect of Temperature Ratio (Tw/T^) on the Aerodynamic Flow Fields for Trajectory Reynolds Number Conditions 

Re^ ia = 0. 22 x 10 6 (Trajectory) ~ D = 257" 

FIGURE 6. SCHLIEREN PICTURES OF THE MODEL AT T /T = 0. 88 AND 0. 22 

w o 


its maximum. A similar increase of is obtained 

by a moderate increase of the Reynolds number. 

Therefore, it is believed that the observed reaction 

of the flow to T /T is caused by effects of the tem- 
w o 

perature ratio on boundary layer transition. The 

results indicate that under proper circumstances 

T /T must be observed as a similarity parameter, 
w o 

VI. VARIABLE POROSITY WALLS 


The MSFC 14 x 14-inch (35. 56 x 35. 56-cm) 
Trisonic Wind Tunnel uses a transonic test section 
with perforated walls. Use of perforated walls to 
reduce the wall interference is an accepted practice. 


Figure 8 shows a typical test section of this type. 
Normally, one fixed type of wall perforation is used 
throughout the transonic range, rather than per- 
forations adapted individually to each Mach number. 
The price for this simplification is a reduced data 
accuracy. Our Facilities Branch has, therefore, 
developed walls of variable porosity (Fig. 9). They 
permit a change of the porosity by simply moving one 
of the two wall plates. The result is illustrated by 
the next two graphs. Figure 10 shows a cone-cylinder 
pressure distribution with fixed-porosity walls set 
for good average cancellation at 1. 0 < M ^ 1. 3; 

Figure 11 shows the measured pressure distribution 
on the same body with our variable porosity walls set 
for optimum cancellation at this Mach number. 

Outside agencies are becoming interested, i. e. , 
LeRC, AEDC. 
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Variation of Drag Coefficient with Temperature 

Ratio (T / T j) f° r Trajectory Reynolds Number Conditions 


FIGURE 7. THE EFFECT OF T /T ON THE DRAG COEFFICIENT 

w o 


VI I. DYNAMIC BALANCE FOR SATURN 
FOREBODY 

The aerodynamically unfavorable nose shape of 
the Saturn vehicles may give rise to dynamic bending 
instabilities. To perform wind-tunnel tests on these 
effects, a sting-mounted dynamic balance (Fig. 12) 
was designed. The balance forms an integral part 
of the model since its frame supports and stiffens 


the light external shell. This balance permits meas- 
urements at preselected frequencies ^ 100 Hz, angles 
of attack, and amplititudes of ± 2 1- degrees about the cen- 
ter of rotation. The frequencies are limited by stresses 
in the escape tower of the model. Preliminary tests 
in the transonic range proved that low internal damping 
was achieved. The aerodynamic damping is 70 to 85 
percent of the total damping measured. The balance 
requires only simple electronic equipment, and permits 
high acceleration due to its ample torque and power. 
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FIGURE 8. PHOTOGRAPH OF 20° CONE- 
CYLINDER IN TUNNEL 


FLOW t=£> 



OPEN CONDITION 


FLOW C=£> 


PARTIALLY CLOSED 
CONDITION 


MOVABLE WALL 


MOVABLE WALL 



FIGURE 9. SKETCH OF VARIABLE POROSITY 
CONCEPT 


FIGURE 12. STING, MOUNTED DYNAMIC BALANCE 
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By 

M. F. Platzer 


SUMMARY 

This study briefly describes some of the problem 
areas in the Saturn-launch-vehicle-development pro- 
gram that are related to unsteady phenomena. The 
two principal types of unsteady phenomena are self- 
excited phenomena and externally -excited phenomena. 
Self-excited phenomena consist of panel flutter, body 
flutter, and fin flutter. Studies conducted in this area 
show that strong influences are exerted by the bound- 
ary layer. Concurrently, studies are being made to 
calculate the nonlinear and the linearized-invisc id- 
flow over pulsating and stationary wavy-walled sur- 
faces. 

The following three methods have been initiated 
for predicting aerodynamic damping: a linearized- 
unsteady-potential equation based upon a generalized 
Adams-Sears iteration process; a linearized method 
of characteristics solution for slowly oscillating 
pointed-bodies of revolution that are of arbitrary 
profile and that are in supersonic flow; and a non- 
linear method of characteristics solution. 

Externally-excited phenomena consist of buffeting 
loads, noise problems, and meteqrological environ- 
ment, such as wind gusts and atmospheric turbulence. 
Research to determine the power-spectral character- 
istics of various buffeting regions is in process. 
Theoretical studies are also being made on the rel- 
ative contribution of regions of high-entropy pro- 
duction in the exhaust stream to the overall intensity, 
propagation, and directivity characteristics of the 
radiated acoustic energy. 

I. INTRODUCTION 

Unsteady aerodynamics are airflows that are 
time dependent and that may produce a dynamic re- 
sponse of the vehicle structure or control system. A 
characteristic feature of these response problems is 
that comparatively small aerodynamic forces may 
produce a very large or even catastrophic dynamic 
response. Aircraft wing flutter and the Tacoma 
Narrows Bridge failure are typical examples. 


Unsteady phenomena can be divided into the general 
classes of self-excited and externally-excited phe- 
nomena. Self-excited phenomena shall consist of: 

(a) Panel flutter. 

(b) Body flutter. 

(c) Fin flutter. 

Externally -excited phenomena shall consist of: 

(a) Buffeting loads. 

(b) Noise problems. 

(c) Meteorological environment, such as wind 
gusts and atmospheric turbulence. 

The primary objective of this study is to describe 
the efforts to provide an accurate description of un- 
steady aerodyanmic forces to be used to furnish input 
information to the structures’ and control systems’ 
designer. No established procedures exist for 
determining an unsteady aerodynamic environment in 
connection with space vehicles or missiles. There- 
fore, this study is primarily of a research nature. 
However, these research results should be converted 
into useful design information as soon as possible. 


II. PANEL FLUTTER 


Panel flutter is the self-excited oscillation of a 
panel in a flow under the action of aerodynamic forces 
generated by the panel motion. The possibility of 
panel flutter was suspected in the early phases of the 
V-2 development. Since then, experiments have shown 
conclusively that flutter of surface skin panels can 
exist. Wind-tunnel tests showed that a traveling, 
wave-type, flutter condition caused panel failure with- 
in a period of a few seconds [ 1] . Recently, panel 
flutter occurred on the X-15 airplane as shown in 
Figures i and 2. In general, structural nonlinearities 
tend to limit the flutter amplitudes and, therefore, 
cause the modes of structural failure to be related to 
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fatigue rather than to the explosive fracture of the 
skin surface. However, it is quite dangerous to 
regard panel flutter only as a fatigue problem. Tests 
of the corrugation stiffened X-15 panels, e.g. , re- 
vealed that flutter was very sudden and very severe. 
The panels failed unless flutter was stopped within 
a period of a few seconds. 



FIGURE 1. REGIONS OF X-15 RESEARCH 
AIRCRAFT AFFECTED BY PANEL FLUTTER 
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FIGURE 2. VARIATION OF PANEL-RESPONSE 
ENVELOPE WITH DYNAMIC PRESSURE OBTAINED 
FROM FLIGHT MEASUREMENTS ON AN X-15 
SIDE -FAIRING PANEL 

For this reason, it was considered necessary by 
the Aero-Astrodynamics and the Propulsion and Ve- 
hicle Engineering Laboratories to test the hat sections 
of the Saturn-IC stage. These tests, which were 
carried out by Boeing personnel in the Ames Research 
Center (ARC) 2- and 11-foot (0. 61- and 3. 35-m) 
tunnels, demonstrated the strong influence exerted 
by the boundary layer. 

Theoretical-flutter calculations are based upon 
inviscid-linearized-flow theory. Only recently has 
the need for taking into account the effect of the 
viscous boundary layer been recognized. For this 
reason, the Unsteady Aerodynamics Branch has 
initiated a cooperative research program with ARC 



and Georgia Institute of Technology to study the in- 
fluence of the boundary layer on panel flutter. A 
decision was made to first explore the stationary flow 
over wavy-walled surfaces for varying boundary layer 
thicknesses. The first test in this program was 
scheduled for January 1966 at the ARC 2-foot (0. 61-m) 
tunnel. The Mach number range to have been explored 
was 0. 8 < M < 1.6. Concurrently, theoretical -inhouse 
studies are being carried out to calculate nonlinear - 
and linearized-inviscid flow over pulsating and sta- 
tionary wavy-walled surfaces [2] . 

III. AERODYNAMIC-DAMPING STUDIES 


As a result of the increased flexibility of the 
launch vehicles in the Saturn class, it has become 
necessary to investigate the aerodynamic damping 
characteristics of these vehicles during flight in the 
first three bending modes. Theoretical approaches 
to predict the aerodynamic damping are especially 
important because the simulation of elastic vehicle 
response in a wind-tunnel test is very complicated, 
and the results of this expensive test may be of limited 
value due to the wind-tunnel noise and the model- 
mount interference. Another factor in favor of the 
theoretical approach is that the elastic model usually 
necessitates a lead time of 1 year. Consequently, 
the test results become available when the vehicle 
design has been frozen. 

For these reasons, the Unsteady Aerodynamics 
Branch tried to formulate theoretical approaches to 
predict the aerodynamic damping. These theoretical 
approaches had to be based upon inviscid -attached - 
flow concepts in order to make the problem mathe- 
matically tractable. So far, the following results 
have been obtained: 

(a) A solution of the linearized unsteady potential 
equation has been worked out based upon a 
generalized Adams -Sears iteration process. 
This work is now complete [ 3] and Figure 3 
shows one result that demonstrates the in- 
fluence of body shape and Mach number upon 
the pitch damping coefficient. 

(b) A linearized method of characteristics 
solution is presently being investigated for 
slowly oscillating pointed bodies of revolution 
that are of arbitrary meridian profile and 
are in supersonic flow. This method is an 
extension of the Oswatitsch-Erdmann char- 
acteristics method. Preliminary results for 
the steady case show a wide applicability 
[4]- 
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(c) A nonlinear method of characteristics 
solution is being investigated by the 
Massachusetts Institute of Technology (MIT) . 
This approach is based upon the superposition 
of small time -dependent angle of attack per- 
turbations on a nonlinear -axisym metric - 
flow field that is determined by the method of 
characteristics. Good results have been 
obtained for the oscillating cone [ 5] . 



FIGURE 3. EFFECT OF MACH NUMBER ON FIXED 
AXIS DAMPING IN PITCH MOMENT COEFFICIENT 
FOR CONVEX AND CONCAVE PARABOLIC OGIVE 
OF THICKNESS RATIO 0. 1 


Calculation of damping characteristics, which 
are based upon attached flow concepts, are not suf- 
ficient for most of the current, launch-vehicle con- 
figurations. Blunt nose cones and steep interstage 
flares, which have become characteristic of space 
boosters, produce large regions of flow separation 
that dominate the aerodynamic loading over the ve- 
hicle in the transonic and the supersonic speed range, 
as shown in Figure 4. As a result, it has become 
necessary to take into account separated flow effects 
upon aerodynamic damping. The success of an 
analytical and experimental method of dynamic anal- 
ysis of the Atlas-Able launch vehicles [6] suggested 
a possible adaptation to the Saturn-launch-vehicle 
configuration. This method is based upon the as- 
sumption that the total lift force on any part of the 
vehicle consists of the following components: 


(a) Local lift, which is lift due to local motion 
at time t. 

(b) Induced lift, which is lift induced by the 
wake created at the nose at time t-At. 


This suggests that local lift is always in phase 
with the lateral vehicle motion and that induced lift 
is out of phase by the time lag At. These two as- 
sumptions make it possible to determine aerodynamic 
damping from static-wind-tunnel tests thus eliminating 
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expensive elastic model tests. This procedure is 
being applied to the various Saturn-launch-vehicle 
configurations. 



FIGURE 4. FLOW FIELD ABOUT TYPICAL 
SATURN SERIES VEHICLE 

IV. BUFFETING 

Buffeting is the pressure oscillation caused by 
separated flow. The launch- vehicle-inter stage areas 
may cause intensive fluctuating pressures, and, 
consequently, severe structural loads as shown in 
Figure 5. In the separated flow region shown in 
Figure 5, a slowly circulating-reverse flow exists, 
overlain by a free-shear flow. The highest fluctuating 
pressure levels are found to occur at the point of 
flow separation and reattachment. Current research 
is devoted to the determination of the power -spectral 
characteristics of these buffeting regions. The Un- 
steady Aerodynamics Branch has several research 
programs to analyze this problem. Perhaps the most 
significant program in this field is the protuberance - 
wind-tunnel test which is being carried out in co- 
operation with ARC. This program is currently in 
the data reduction process and the results should be 
available in the near future. 



FIGURE 5. TYPICAL FLUCTUATING PRESSURE 
ENVIRONMENT FOR VEHICLE INTERSTAGE TYPE 
GEOMETRIES 
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V. ENGINE-GENERATED NOISE 


Research in the field of engine-generated acous- 
tics has been directed towards the determination of 
the predominant noise generation mechanisms in 
high- temperature, rocket-exhaust flows. Present 
investigations include theoretical studies of the 
relative contribution of regions of high entropy pro- 
duction in the exhaust stream to the overall intensity, 
propagation, and directivity characteristics of the 
radiated energy. This work is being supplemented 
by experimental rocket exhaust studies using a 
substitute gas to simulate the predominant flow 
characteristics of the H-l engine at a considerably 
lower temperature, thus allowing measurements to 
be made in the flow. 

VI. GROUND-WINDS PROBLEM 


A launch vehicle erected on the launch pad be- 
fore lift-off is exposed to surface winds that may 
generate large dynamic loads under certain circum- 
stances. Downstream from the body, a wake flow is 
created that is Reynolds number dependent and has 
characteristics which are difficult to describe math- 
ematically. In certain cases, the well known Karman 
vortex street is obtained. At the higher Reynolds 
numbers, the vortex shedding is very dependent upon 
the type of flow in the vicinity of the body ( laminar - 
or turbulent-boundary layer) . Therefore, a dis- 
tinction must be made between subcritical and super- 
critical flow regimes. Insufficient understanding of 
the mechanisms involved and lack of reliable cal- 
culation methods force most investigators to extract 
the aerodynamic forcing function from an analysis of 
the aeroelastic response of various launch -vehicle 
configurations. The Unsteady Aerodynamics Branch 
has initiated an extensive wind-tunnel program on 
the different Saturn I, IB, and V configurations. Most 
of these tests are conducted at the Langley Research 
Center (LRC) and require aeroelastic models of 
about 5-percent scale. Unfortunately, for the Saturn 
V configuration, a full-scale Reynolds number cannot 
be simulated, thus introducing a further degree of 
uncertainty. Typical vehicle responses are shown 
in Figures 6 and 7. In both figures the vehicle base 
bending moment, due to steady drag, dynamic drag, 
and dynamic lift, is plotted versus the wind velocity. 

A typical noncritical response is shown in Figure 6, 
while a typical critical response due to dynamic lift 
at a specific wind velocity is shown in Figure 7. 


The critical response is of primary concern to 
the vehicle designer since it may impose exceedingly 
high loads upon the vehicle. Unfortunately, the reason 
for the occurrence of such a critical response is 
little understood. The possibility for self-excited 
vibrations cannot be excluded at this time. There- 
fore, a cooperative research program between LRC, 
MSFC, and the Martin Company was formulated to 
investigate the aerodynamic forces on oscillating 
two-dimensional cylinders. These tests were con- 
ducted recently at LRC and the data are in the process 



FIGURE 6. TYPICAL NON-CRITICAL RESPONSE 



FIGURE 7. TYPICAL CRITICAL RESPONSE 
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of evaluation. This experimental program is being 
supplemented by an analytical studies program that 
should provide a better understanding of the basic 
flow-field characteristics. 


Unsteady aerodynamic effects present the launch- 
vehicle designer with a series of important and 
challenging problems. Since these phenomena are 
little understood, a continuous and vigorous research 
effort is required. 


VII. CONCLUSIONS 
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BASE-HEATING RESEARCH REVIEW 


By 

Homer Wilson 


SUMMARY 


The purpose of the studies described in this re- 
port is to gain more information related to the causes 
and control of heating problems. To help gain this 
information, an external-flow facility was built to 
enable the Reynolds number to be varied by almost 
an order of magnitude; and at the same time, to 
provide a facility that would eliminate costly large- 
scale, wind-tunnel tests. Data obtained from the 
external-flow facility compared reasonably well with 
large-scale, wind-tunnel tests. By using this facility, 
some limited Reynolds numbers on base heating were 
determined. A similar, but larger and more ver- 
satile, facility has been proposed by MSFC. 

This study also indicated that a flow field pre- 
vented the radiation of vehicles with clustered engine 
configurations from decreasing to near zero radiation 
as in single engine vehicles. The flow field results 
from an increasing altitude and decreasing ambient 
pressure causing the exhaust jets from the engines 
in a cluster configuration to expand and to impinge 
with the other jets in the cluster forming a shock 
wave. The shock waves force the exhaust gases 
back toward the base plate of the engine. The reverse 
exhaust gases will then flow radially outward between 
the nozzles into the atmosphere and may choke the 
engine at sufficiently high altitudes. This condition 
has caused severe heating on the Polaris and the 
center star areas of the Saturn vehicles. 

For clustered LOX-H 2 upper-stage configurations, 
radiation may not be the predominant cause of heating; 
but it should be accurately computed since the stages 
fly a long period of time. In order to compute the 
radiation, the ’’band model” concept was used. To 
use this concept, it is necessary to know parameters 
such as spectral-line strength, spectral-line spacing 
and line half-width. Within the next year, all needed 
parameters are expected to be available. 

Another problem studied was the dumping of 
combustibles along the vehicle or near the base area. 

If fuel-rich combustibles are entrained in the base 
region, large heating rates may occur. The results 
from the study show that upstream ejection of 


hydrogen can be accomplished without hazard of 
ignition or combustion by venting parallel to the 
surface through a perfectly expanded nozzle if it is 
displaced far enough from the vehicle to reduce hy- 
drogen levels below flammability limits. 


I. INTRODUCTION 


The purpose of a research program in heating is 
to gain information which will improve the design of 
future space vehicles. The base heating problem on 
rocket vehicles is a combination of heating by con- 
vection and by radiation. For LOX-RP powered 
launch vehicles, radiation is high at sea level and 
decreases at higher altitudes (Fig. 1). The radiation 
at higher altitudes for single -engine vehicles like 
Jupiter decreases to near zero. Although the radiation 



Choked-Flow Condition 


120 

1000ft 


Altitude, 


Engine Vehicle 


Vehicle 


Altitude, 1000 ft 


FIGURE 1. VARIATION OF RADIATIVE AND CON- 
VECTIVE BASE HEATING WITH ALTITUDE FOR 
VARIOUS CONFIGURATIONS 
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of vehicles with clustered engine configurations is 
less at higher altitudes than at sea level, the im- 
pinging regions at higher altitudes apparently bring 
about a significant increase in the radioactive flux 
keeping the radiation of the vehicle from nearing 
zero. Upper stages powered by LOX-H 2 engines 
have little radiation, but the radiation should be ac- 
curately computed since the vehicles fly over a long 
period of time. 

Convective heating for single -engine vehicles 
generally is small at all altitudes. For vehicles 
with clustered engine configurations, the impinging 
of the jets at high altitude sets up a recirculating 
flow field which causes large heating within the en- 
closed area. 

Another problem to be reckoned with is the 
dumping of combustibles along the vehicle or near 
the base area. If fuel-rich combustibles are en- 
trained in the base region, large heating rates may 
occur. These are usually a maximum in the tran- 
sonic region and disappear above Mach 3 due to the 
absence of oxygen to support combustion. 


II. SHORT-DURATION TECHNIQUE 


Over the years, a wealth of flight and model 
data has been accumulated which was used as esti- 
mates of the thermal criteria for design purposes. 
Because these estimates have not been satisfactorily 
substantial, time has been spent on programs aimed 
at improving methods of generating design data and 
getting a better understanding of the problems. During 
early Saturn days, model data were obtained from 
scaled, long-duration, model tests. This type of 
testing proved costly, complex, and generally pro- 
vided data over a limited range of variables. About 
4 years ago, the possibility of using other techniques 
for obtaining model data was explored and the short- 
duration technique was developed. This technique is 
illustrated in Figure 2. 

Using this technique, the propellants are stored 
in a pair of charge tubes and are contained by a pair 
of diaphragms or quick-acting valves at the end of 
each tube. When the diaphragms are ruptured, the 
gases flow into an injector, then into the combustion 
chamber where burning takes place. Ignition is 
usually by a spark plug. For high-altitude testing, 
the model is mounted in a vacuum tank. Testing 
time is determined by the time required for an ex- 
pansion wave to traverse the charge tube (charge- 
tube length) , the time required for the blast wave to 
travel to the tank walls and return, or the heating 
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limitations of the model. Early results on a base 
plate indicate that thermal and pressure fields were 
established in \ to 1 millisecond and that several 
milliseconds of steady flow are available which is 
ample time to obtain the data. 



FIGURE 2. COMBUSTOR SCHEMATIC, SUPPLY 
TUBE-TANK WAVE DIAGRAM AND INSTALLATION 

The short-duration technique has been used to 
obtain design data for S-IV, S-II, S-V, S-I, and 
Centaur stages. On a firing- to -firing basis, the 
short-duration technique is cheaper by an order of 
magnitude than the long-duration technique and pro- 
vides data that compare reasonably well with long- 
duration model and hi gh- altitude , flight data. The 
big advantage of the short-duration technique is that 
many variables can be easily studied. The effects of 
engine deflection on base heating is shown in Figure 
3a for a typical, five-engine configuration. The re- 
sults show that on deflecting two outboard engines 
toward each other by 4 degrees the maximum heating 
rate increases by a factor of two. This is particularly 
important for design of the upper stages such as S-U. 

Some effects of nozzle-wall temperature on base- 
heating rate is shown in Figure 3b. Increasing the 
nozzle-wall temperature from 80 to 1000° F (26.67 
to 537. 78° C) results in a significant increase in 
base-heating rate. This results from the fact that 
the low-energy, exhaust gases that are reversed in- 
to the base have passed along the cooled nozzle-wall 
temperature. 

The effect of a temperature step or discontinuity 
in the measurement of the film coefficient is shown 
in Figure 4a. Cooling the gage results in a large 
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FIGURE 3. DATA OBTAINED USING SHORT 
DURATION TECHNIQUES OF VARIOUS ENGINE 
CONFIGURATIONS 

increase in the film coefficient (almost a factor of 
two) as compared to a case where the gage tem- 
perature is equal to the plate temperature. This 
indicates that large errors may be present in our 
flight measurements where the copper-slug calori- 
meter temperature may be significantly different 
from that of the heat shield. 

The effect of pressure of Reynolds number on 
base-film coefficient and recovery temperature is 
shown in Figure 4b. The initial portion of the curve 
shows that the film coefficient varies proportionally 
to the \ power of the pressure as would be expected 
for laminar flow, and the latter part of the curve 
shows that the film coefficient varies proportionally 
to the 4/5 power of the pressure as would be expected 
for turbulent flow. The base-recovery temperature 
varies proportionally to the 1/5 power of the pressure. 
The dependency of the recovery temperature on pres- 
sure or Reynolds number is brought about by the 
fact that the gases that get reversed into the base area 
are contained within the boundary layer of the nozzle, 
the state or energy level of which is primarily de- 
termined by the Reynolds number. 

1 1 1. EXTERNAL-FLOW REYNOLDS NUMBER 
PROGRAM 


In comparing model to flight data, where ambient 
flow effects are important, several discrepancies 
have occurred which seem to indicate that Reynolds 
number effects are present. Approximately 2 years 
ago, a program was initiated to build an external- 
flow facility which would enable the Reynolds number 
to be varied by almost an order of magnitude and at 
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the same time provide a facility that would eliminate 
costly large-scale, wind-tunnel tests. The facility 
has been developed and has been in operation for 
about 1 year (Fig. 5) . It consists of a tube approx- 
imately 4 feet (1. 22 m) in diameter and 30 feet 
(9. 14 m) in length in which high pressure gas, 200 
psi ( 1378951. 44N/m 2 ) and lower is stored. A dia- 
phragm retains the gas at the end of the tube. To 
the right of the diaphragm is the convergent-diver- 
gent nozzle which allows the Mach number to be varied. 
The flow is exhausted into a vacuum tank which is 
located on the far right of Figure 5. 



FIGURE 5. SCHEMATIC OF EXTERNAL FLOW 
FACILITY WITH SOME TYPICAL PRESSURE TRACES 

To perform a test, the diaphragm is ruptured 
and a flow of a certain pressure level and Mach 
number in the nozzle is obtained for a few milli- 
seconds. About the same time that steady flow is 
established in the nozzle, the model is fired and the 
test takes place. Hot and cold flow tests have been 
made in this facility. Figure 5 also shows typical 
traces from the early tests at a Mach number of 2. 0. 

P g is the shroud pressure just ahead of the base; P^ 

is the stagnation or tube pressure just downstream of 
the diaphragm; P^ is the combustion chamber pressure 

within the model; P^ is the typical base pressure 
trace; and P^ is the pressure within the vacuum tank. 
The large variation in P^ at approximately 50 milli- 
seconds indicates the return of the blast wave and 
the end of the test. 

The results show that for this length of storage 
tube approximately 20 milliseconds of constant-flow 
conditions are available for making tests. Due to the 


late firing of the combustor, only about 10 milli- 
seconds were used. The storage tube and combustor 
supply tubes have been lengthened so that considerably 
longer testing time can be obtained. 

Data obtained from the external-flow facility com- 
pared reasonably well with large-scale, wind-tunnel 
tests. The external-flow facility has been built and 
put into operation for a cost of less than $75, 000. By 
using this facility, some limited Reynolds number 
effects on base heating were determined. This 
facility may reduce the number of, or even eliminate, 
costly wind-tunnel tests for Saturn base-heating 
studies. A similar, but larger and more versatile 
facility, has been proposed by MSFC and some fea- 
sibility tests for the facility have been made by MSFC. 

IV. FLOW-FIELD VISUALIZATION 


The flow field results from an increasing altitude 
and decreasing ambient pressure causing the exhaust 
jets from the rocket engines in a cluster configuration 
to expand and to impinge with the other jets in the 
cluster forming a shock wave. Figure 6a shows a 
sketch of the flow field. Because of the formation of 
the shock waves and the resulting pressure rise, a 
portion of the low-energy, exhaust gases within the 
boundary layer of the nozzle wall and in the mixing 
layer around the periphery of the jets will not have 
sufficient energy to negotiate the pressure rise and 
will be forced back toward the base plate of the engine. 
The reverse exhaust gases will then flow radially 
outward between the nozzles into the atmosphere, or 
ambient flow, and may choke the engines at sufficiently 
high altitudes. This condition has caused severe 
heating on the Polaris and the center star areas of the 
Saturn vehicles. 

Figure 6b is a schlieren photograph of the shock 
waves formed by the impinging jets. Some results 
from electron-beam, density measurements on the 
short-duration, S-IV four -engine model are shown in 
Figure 6c. The solid lines represent lines of constant 
density in the form of density ratio ( Py^P Q ) where 

p^ is the ambient density in the reverse flow field 

and is the density within the combustion chamber. 

The gas is hot nitrogen (2500° R) at a pressure of 
200 psi ( 1378951. 44N/m 2 ) . The dashed lines repre- 
sented the nozzle and the approximate plume impinge- 
ment. The abscissa is the base radius and the ordinate 
is the number of inches downstream of the base plate. 
Starting at the centerline and proceeding toward the 
base plate, the density decreases which indicates an 
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SKETCH OF FLOW FIELD 


SCHLIEREN PHOTOGRAPH 
OF SHOCK WAVE 


DENSITY FIELD BASED ON ELECTRON 
BEAM MEASUREMENTS 
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FIGURE 6. RESULTS FROM FLOW- FIELD, VISUALIZATION STUDY 


accelerating flow field. A short distance above the 
base plate, a large increase in density occurs which 
indicates the presence of a shock wave. Moving 
farther from the base plate, the density decreases 
indicating an accelerating flow in that direction. This 
type of flow field has been studied and verified with 
conventional-pressure instrumentation. 

Only some of the initial results from the electron- 
beam tests are shown in Figure 6. Hardware has 
been assembled and an attempt at measuring tem- 
perature with the electron beam will be made in the 
near future. 

An instrument for the measurement of local gas 
velocity using the optical Doppler technique is also 
under development. The Doppler frequency shift 
depends on the velocity of the moving particles and 
on the geometry of the scattering media. If the 
scattering geometry is fixed, measurement of the 
Doppler shift gives sufficient information to determine 
the velocity of the moving particles. This instru- 
ment has been built and velocity measurements have 
been made up to 150 feet per second (45. 72 m/s). 

The instrument is presently being installed in a wind 
tunnel for measurements of Mach 3 to 5. This 


technique seems applicable also for gas -turbulence 
measurements , dynamic-pressure-transducer- 
calibration, panel flutter studies, and many others. 

V. RADIATION 

As stated earlier, radiation is the main con- 
tribution to base heating for a properly designed LOX- 
RP first stage. For clustered LOX-H 2 upper-stage 
configurations, radiation may not be the predominant 
cause of heating; but it should be accurately computed 
since the stages fly a long period of time. In order 
to compute the radiation, the flow-field (velocity, 
pressure, temperature, composition, etc.), ab- 
sorption data for radiating components and a radiant 
program for calculating the intensities must be known. 
Programs are in existence for calculation of the flow 
field from a single-engine configuration for frozen 
or equilibrium flow. Some rough estimates of the 
flow fields in the impinging regions are currently 
being used for design purposes. Recently, two pro- 
grams have been started for computing the flow field 
for clustered configurations. One is a three-dimensional 
method of characteristics at Norair and the other is a 
finite difference scheme at General Dynamic s/Convair. 
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Until some new programs were initiated, the 
techniques and data available to calculate the radiant 
heat from H 2 O f C0 2 , and CO used the total emmissivity 
data of Hottel measurements during the late twenties 
and early thirties. These data could not be accurately 
extended to the range of low pressure, high tem- 
peratures, and nonisothermal conditions existing in 
the fuels of current interest, RL-10, J-2, etc. It 
is possible to compute some of the data needed with 
existing high-speed computers, but this is difficult 
and would require an enormous amount of computing 
time. Therefore, the "band model” concept was used 
to describe the absorption coefficient for the radiating 
components. To use this concept, it is necessary to 
know such parameters as spectral- line strength, 
spectral-line spacing, and line half- width. Programs 
have been initiated and some results have been ob- 
tained for some of the band-model parameters needed. 
Within the next year, it is expected that all para- 
meters needed will be measured so that radiation 
from present motors can be computable within the 
accuracy of existing flow fields. 

Figure 7 presents some early results of the 
"band model” concept. The upper part of the curve 
compares experimental data with "band model" data 
for water vapor. The results indicate that the "band 
model” concept is relatively good. The graph in the 
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FIGURE 7. COMPARISON OF EXPERIMENTAL 
DATA WITH "BAND MODEL” RESULTS FOR WATER 
VAPOR 


lower left-hand corner compares the recent emmission 
data of water vapor of GD/C (Ferriso and Ludwig) 
with Hottel's results. The results indicate that de- 
sign data based on Hottel T s results would under- 
estimate heating rates. These results are supported 
by early experimental data. A comparison of the 
experimental absorption data for carbon with theory 
(Stull and PI ass, modified for index of refraction) is 
shown in the lower right-hand graph. The results of 
the experimental absorption data are relatively good 
except at the lower wavelengths. 


VI. INVESTIGATION OF COMBUSTION OF 
HYDROGEN IN A HYPERSONIC AIRSTREAM 


Heat leaks to hydrogen- filled fuel tanks made it 
necessary in some cases to dump combustible hy- 
drogen during flight. Since hydrogen will mix with 
the ambient air flow, the possibility of it combusting 
and releasing large amounts of heat in the vicinity of 
the vehicle has been investigated through a contract 
with General Applied Sciences Laboratories. Chemical 
kinetic studies were performed to determine where 
in the launch trajectory conditions would be favorable 
for combustion. Fluid mechanical mixing models 
were set up to investigate two dumping configurations: 
tangential slot injection, and injection by means of 
jets located parallel to, but removed from, the sur- 
face of the vehicle (Fig. 8). Also more sophisticated 
flow models employing finite-rate chemistry kinetics 
were developed. 



FIGURE 8. MIXING OF HYDROGEN JET WITH 
VEHICLE BOUNDARY LAYER 

To help determine where in the launch trajectory 
conditions would be most favorable for combustion, 
two-phase flow phenomenon and the condensation and 
evaporation associated with the dumping of cryogenetic 
hydrogen were investigated. The analytical formu- 
lations from the above investigation were also sup- 
ported by experimental investigations of both tan- 
gential slot and jet' dumping configurations. The re- 
sults from the above investigations have shown that 
the hydrogen and air will mix in lengths which are 
short compared to practical vehicle dimensions and 
will form mixtures whose composition will fall within 
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the combustion limits of the air-hydrogen system. 

If local hot spots are present in the flow field, such 
as high-static temperature occurring in the boundary 
layer or increased heating due to protuberances, 
combustion can be initiated. Combustion may not 
begin immediately, but may be suppressed because 
of the high heat capacity of the hydrogen. 

The results of the investigation also show that 


upstream ejection of hydrogen can be accomplished 
without hazard of ignition or combustion by venting 
parallel to the surface through a perfectly expanded 
nozzle if it is displaced far enough from the vehicle 
to reduce hydrogen levels below flammability limits 


Some results of this program are being used in 
the Saturn Centaur studies. 
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RAREFIED-GAS DYNAMICS 


By 

James O. Ballance 


SUMMARY 


This study develops two new techniques for the 
analysis of complicated impact pressure probes at 
various speed ratios and angles of attack. The first 
technique is based on the analogy of radiant-heat 
transfer to that of free-molecular-flow. The second 
technique uses the Monte Carlo method of selecting 
random numbers to follow the paths of individual 
molecules in a system. 

The calculation of the drag coefficient of a con- 
cave hemispherical shell is extended from the 
limiting case where all incoming flow is perpendicular 
to the opening plane of the shell to the case where the 
incoming flow is at some angle a to the opening plane 
of the shell. 

A modification of the Monte Carlo method is 
being used in the flow analysis of a cylindrical duct 
in the transition-flow regime. This method il- 
lustrates that a Knudsen number value based only on 
the diameter of the duct is not sufficient to define the 
limit for free -molecular -flow but that the length of 
the duct also has an important effect on the flow. 

Thermal-accommodation coefficients have only 
been measured for low-energy collisions and even 
then the results are questionable. A program is being 
conducted to measure the momentum and energy ac- 
commodation coefficients up to and exceeding the 
escape velocities on typical engineering surfaces. 

The experimental facility used to conduct these 
measurements is nearly completed, and it is hoped 
that useful data will be obtained by the end of the 
year. 


I. INTRODUCTION 


In the discussion of this topic, it is necessary to 
establish a reference. The best reference is the mean 
distance that a gas molecule travels between collisions 


with other gas molecules (mean-free path). For 
standard conditions, the mean-free path of an air mole- 
cule is approximately 6. 6 x 10" 8 meters. As the num- 
ber of molecules per cubic meter (number density de- 
creases, the me an- free path increases. Figure 1 
shows the mean-free path for several different values 
of number density. To obtain a better understanding 
of the physical significance of these values of mean- 
free path, the altitudes for the various density levels 
as given by the United States Standard Atmosphere 
( 1962) are also shown in Figure 1. When the number 
density is very low, the mean-free path is so large 
that the collisions between gas molecules can be 
ignored for most studies involving conventional sized 
objects. This region is known as the free molecular flow 
regime. As the number density increases, the mean- 
free path becomes very small, molecular collisions 
increase, and the flow approaches the continuum- 
flow regime. The parameter which describes these 
regions with reference to the object under study is 
the Knudsen number (Kn) . The Knudsen number is 
the ratio of the mean-free path to the characteristic 
dimension of the object. For example, the charac- 
teristic dimension of a cylinder is normally taken to 
be the diameter of the cylinder. Figure 2 shows the 
flow regimes as normally defined in terms of the 
Knudsen number. 


(US 

Standard 

Atmosphere , 

1962 ) 



ALTITUDE 

NUMBER 

DENSITY 

MEAN 

FREE PATH 

(kilometer) 

(molecule 

/ m 3 ) 


(meter) 

0 

2.55 

X 

10+25 

6.6 

X 

10 - 8 

50 

2.14 

X 

10 +22 

7.9 

X 

10 - 5 

100 

1.04 

X 

10 + 19 

1.6 

X 

10 - 1 

150 

4.10 

X 

10 + 16 

4.0 

X 

10 + 1 

200 

7.82 

X 

10 + 15 

2.2 

X 

10 + 2 

250 

2.49 

X 

10+15 

6.8 

X 

10 + 2 

300 

9.50 

X 

10 + 14 

1.8 

X 

10 + 3 

350 

4.15 

X 

10 + 14 

4.1 

X 

10 + 3 

400 

1.96 

X 

10+14 

8.6 

X 

10 + 3 


FIGURE 1. MEAN-FREE PATH OF AIR AT 
VARIOUS ALTITUDES 
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KN = KNUDSEN NUMBER 

= X/L WHERE X = MEAN FREE PATH 

L = CHARACTERISTIC DIMENSION 


FREE MOLECULAR FLOW 

KN 

> 10 

TRANSITION FLOW 

10 

> KN > 0.1 

SLIP FLOW 

0.1 

> KN 


The first method uses the analogy of radiant- 
heat transfer and free -molecular -flow. In free- 
molecular-flow, an assumption is made that molecules 
leave a surface element in a direction proportional to 
the cosine of the angle between the direction vector 
and the normal to the surface. Also, the molecules 
are assumed to travel in rectilinear paths between 
collisions with the surfaces. These assumptions are 
quite similar to those used in radiant-heat- transfer 
calculations. Figure 3 compares the molecular flow 
rate through a duct calculated by this method to the 
numerical solutions of Clausing [ 1] . 


FIGURE 2. FLOW REGIMES FOR RAREFIED-GAS- 
FLOW 

This study is mainly concerned with the free- 
molecular-flow (Kn > 10) and the near -free - 
molecular-flow (1 < Kn < 10) portion of the transition- 
flow regime. This study is grouped into the following 
four basic types: 

(a) Pressure probe characteristics for en- 
vironmental measuring devices. 

(b) Aerodynamic coefficients for orbiting ve- 
hicles in both transitional -flow and free- 
molecule-flow regimes. 



(c) Jet spending characteristics for propulsive 
jets, vent systems, etc., at orbital altitudes. 

(d) Surface physics studies for high-velocity, 
gas -surface interaction. 

II. PRESSURE PROBE CHARACTERISTICS 


Measurements of density, composition, tem- 
perature, and other thermodynamic properties of 
the upper atmosphere are usually made on a gas that 
has been brought to rest in a density and composition- 
sensitive gage connected to the free stream by an 
orifice or an orifice-restricted tube. The flow param- 
eters that relate the properties in the gage with the 
free stream are a function of the orifice size, tube 
length to diameter ratio, energy-accommodation 
coefficients, relative velocities of the tube and free- 
stream, thermal motion, angle of attack, etc. 
Calculated solutions for simple, impact-pressure 
probes at various speed ratios and angles of attack 
have been made; but when more complicated geo- 
metries are considered, no convenient solutions 
exist. Two new methods are being used for this 
analysis which not only yield the desired parameters, 
but also give better insight into the problems. 


FIGURE 3. MOLECULAR- FLOW RATE THROUGH 
DUCT AS A FUNCTION OF DUCT LENGTH TO 
RADIUS RATIO, L/R 

A second method uses Monte Carlo techniques 
where random numbers are used to follow the path 
of individual molecules through the system. By 
following a large number of molecules (e. g. , 10, 000) , 
the properties of the entire system can be determined. 
Figure 4 shows the transmission probabilities for 
tubes as a function of speed ratio. The speed ratio 
is the ratio of gage velocity to the random -thermal 
speed of the gas molecules. The data in this study 
are for cases where there is relative motion between 
the tube and the gas, such as a density gage on a 
satellite. These data are believed to be the first ever 
published for orifice -restricted tubes. Figure 5 
shows transmission probabilities as a function of the 
angle of attack. Here the speed ratio (S) is 2 and the 
length- to- radius ratio of the tube is 4. The lower 
curve is for the case where radius of the orifice 
(Ro) is the same as the diameter of the tube (that 
is, there is no restriction) . The upper curve shows 
the case where Ro is 0. 707 of the radius of the tube 
(R). Also shown in Figure 5 are calculated values 
for the same case using the techniques developed in 
"Transmission Probability Determination with 
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Directed Mass Motion and with Mean Free Path Con- 
sideration" [2], 



FIGURE 4. TRANSMISSION PROBABILITY FOR 
TUBES AT ZERO ANGLE OF ATTACK FOR VARIOUS 
SPEED RATIOS 



unusual configuration (the Pegasus satellite) is 
being made. From this study, more realistic drag 
coefficients will be generated. 

The calculation of the drag coefficient of a 
concave -hemispherical shell has been extended from 
the limiting case where all incoming flow is per- 
pendicular to the opening plane of the shell to the 
case where the flow is at some angle a with respect 
to the plane of the shell [ 3) . Figure 6 presents the 
drag coefficient of a concave hemisphere for speed 
ratios of 4, 6, and 10 where the temperature ratio 
of the surface of the hemisphere (T^) to that of the 

incident stream (T) is 0. 4. Zero angle of attack is 
when the incoming flow is perpendicular to the plane 
of the opening of the hemisphere. Using Monte 
Carlo methods, this effort is being duplicated so 
that the assumption of hyperthermal flow may be 
investigated. Hyperthermal flow is that flow which 
exists when the ratio of the relative speed of an ob- 
ject to the random thermal motion maybe ignored, 
e. g. , S = 6. The flow may be considered to be a 
beam of molecules traveling in parallel paths at a 
constant velocity. Present theories extend this 
type of analysis to speed ratios as low as 6. 



FIGURE 5. TRANSMISSION PROBABILITIES FOR 
ORIFICE RESTRICTED TUBE AT VARIOUS 
ANGLES OF ATTACK 

III. AERODYNAMIC COEFFICIENTS 


Aerodynamic coefficients, such as lift or drag, 
are easily calculated for free -molecular-flow on 
simple convex bodies; but when complex or concave 
bodies are considered, there are no satisfactory 
solutions. To illustrate the problem involved and to 
analyze an orbital system, a thorough study of an 


FIGURE 6. DRAG COEFFICIENT, T /T = 0.4 

w 

IV. JET-SPREADING CHARACTERISTICS 


The problem of gas impingement of surfaces 
during low-mass flow or low-density venting has 
necessitated a closer look at the distribution of 
flow from tubes. While there is a great deal of in- 
formation available for the situation where the gas 
(both in the tube and the free stream) is in the free- 
molecule-flow condition, little information is 
available for the case where the gas in the system 
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is in the transition-flow regime. A modification of 
the Monte Carlo method is again being used for this 
study. Figure 7 presents typical data where this 
modified method was used to study the transition- 
flow regime in a cylindrical duct. The probability 



FIGURE 7, NORMALIZED TRANSMISSION 
PROBABILITIES AS A FUNCTION OF KNUDSEN 
NUMBERS FOR I/A VALUES OF 0. 2, 0. 5, 1. 0 
AND 2. 0 

(K ) is the value deduced by Clausing [ 1] for the 
case of an infinitely long, mean-free path. The 
transmission probability (K ) is obtained from the 

A 

modified method for some mean-free path X . When 
the ratio K /K is equal to 1, the system can be as- 

A C 

sumed to be in free-molecular-flow. The four curves 
shown in Figure 7 are for four different length-to- 
radius ratios (L/A). The curves illustrate that a 
Knuds en number value based only on the diameter of 
the duct is not sufficient to define the limit for free- 
molecular-flow but that the length of the duct also 
has an important effect on the flow. 

V. SURFACE-PHYSICS STUDIES 


While the gross effects of studies in rarefied 
flow are fairly well bounded, they are valid only when 
consideration is given to the interactions of the gas 
and the surfaces of the system under study. The 
microscopic actions in these flows are little known 
and understood. The parameters that specify these 
interactions are the reflection and accommodation 
coefficients. The reflection coefficients are: 

T. _ T 

i r 

o = 

T. 

1 



where 

<t = reflection coefficient due to shear stress 

cr ? = reflection coefficient due to pressures 

t = shear stress due to incident mass 
i 

t = shear stress due to reflected mass 
r 

P = pressure due to incoming free-stream-molecules 

]? r = pressure due to reflected molecules 

P = pressure due to reflected molecules at sur- 
face temperature 

When a - <r' = 0 specular reflections occur and when 
<j = cr f = 1 diffuse reflections occur. The thermal 
accommodation coefficient is: 

dE. - dE 
i r 

a = 

dE. - dE 
1 s 

where 

dE^ = energy due to incoming free-stream-molecules 

dE^ = energy due to reflection molecules 

dE = energy due to reflection molecules at surface 
temperature 

When a is equal to zero, there is no exchange of 
energy and when a is equal to one, perfect accom- 
modation takes place. A good method exists for 
measuring the momentum accommodation of gas 
molecules on surfaces under study. Thermal- 
accommodation coefficients have only been measured 
for low-energy collisions and even then the results 
are questionable. A program is being conducted to 
measure the momentum and energy accommodation 
coefficients up to and exceeding the escape velocities 
on typical engineering surfaces. The experimental 
facility is nearly completed, and it is hoped that use- 
ful data will be obtained before the end of the year. 
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VI. EXPERIMENTAL FACILITIES 

To verify the limits for the analytical solutions 
to the problems in this study, it is necessary to 
perform some experiments. These experiments may 
be performed in a low-density-research chamber 
equivalent to that in use by the Aerodynamics Divi- 
sion. This facility is 3. 5 ft. ( 1. 07 m) diameter by 


14 ft. ( 4. 27 n) in length. A nozzle is being designed 
to produce a flow of Mach 4. However, it is doubtful 
that experimental facilities will ever be developed to 
adequately explore all the problems outlined in this 
study. 

The Aerodynamics Division is conducting a 
modest program that promises to make a significant 
contribution in the field of rarefied-gas dynamics. 
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TURBULENT FLUCTUATION MEASUREMENTS 
WITH THE CROSSED-BEAM METHOD 

By 

F. R. Krause, M. J> Fisher,* and R. E. Larson** 


SUMMARY 

A new optical method has been developed for 
studying aspects of launch- vehicle turbulence. The 
heart of the method is a new test arrangement for 
the remote sensing of local changes in radiative 
power. Early experimental results are given and 
show that most statistical turbulence parameters can 
be approximated which are commonly derived from 
n two point product mean values. " The measurement 
of local thermodynamics’ properties in turbulent 
flows becomes conceivable. 

I. INTRODUCTION 


The measurement of turbulent fluctuations in 
the environment of launch vehicles is necessary to 
provide the input for the statistical analysis of re- 
sponses like rigid-body motions, elastic deformations, 
as well as heat transfer and acoustical loads. Solid 
probes introduce stabilizing walls and/or additional 
shock waves which interfere severely with the fluc- 
tuations of interest. A further instrumentation 
problem is to find probes with a sufficiently high 
resolution in space and time so the fluctuations are 
not integrated out. So far, all ground tests and flight 
instrumentation are restricted to pressure transducers, 
thin film gages, and accelerometers that are flush 
mounted in the vehicle surface. Fluctuation measure- 
ments in the flow have not proved to be feasible. 

The required inputs or aerodynamic -forcing 
functions are mathematically similar insofar as they 
should be expressed by an area integral over a space- 
time correlation function. The inherent numerical 
integration requires a large number of points and the 
experimental and numerical effort in providing these 

* IIT Research Institute 

** Applied Science Division, Litton Industries 


pointwise estimates is almost prohibitive, even if a 
large number of transducers are used in the vehicle 
wails. Therefore, this study proposes to use the 
integrating features of optical beams for a ’’one shot” 
estimate of area-integrated-correlation functions in 
the flow. 

1 1. OPTICAL INTEGRATION OVER CROSS- 
CORRELATION AREAS 


Statistical thermodynamics show that the light 
extinction coefficient p is determined by the thermo- 
dynamic properties of the flow [ 1] . Therefore, its 
space -time correlation might be used to study the 
turbulent fluctuations of thermodynamic flow prop- 
erties. 

Consider the experimental arrangement in 
Figure 1. Two narrow beams of white light traverse 
a jet in the y and z directions. Using an optical 
filter in front of the photo detectors, the received 
radiative power is limited to the wavelength interval 
A A centered around the wavelength A. The intensity 
fluctuations of the y and z beams are now correlated 
and normalized with the mean values. These oper- 
ations lead to the measurable quantity 

< i (t) i (t + t )> 

Z V 

G(|, t) = — * — (1) 

y z 

Fisher and Krause [2] show analytically that the 
measurable quantity is related to the two-point, 

product-mean value of the light extinction co- 
efficient 

%U.T) X = (2) 

.<H' (s, y, z, t)fi' (x + j, y + tj, z + £, t + T)> 
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The crotsed beam technique gives the integrated cross correlation function of light extinction coefficients. 

♦© 
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FIGURE 1. OPTICAL INTEGRATION OVER WAVE FRONTS IN TURBULENT FLOWS 


by an area integration across the plane in which the 
beams have been aligned. 

+ oo 

G <«’ T > X = / / « r) x dr} d? . (3) 

— OO 

The derivation of equation (3) shows that the wanted 
"one shot” estimate of area integrated correlation 
functions can be obtained, even in statistically in- 
homogeneous media, as long as there is one direction 
of homogeneity Left. The wavelength A and the spec- 
troscopic resolution AA may be adjusted to a specific 
thermodynamic property, such as density, species 
concentrations, and, hopefully, temperatures. 

III. APPROXIMATION OF POINT 
MEASUREMENTS 


The random nature of turbulent fluctuations as- 
sures that the integrand in equation ( 3) drops to zero 


over finite and mostly small distances. In this 
fashion, local information can be obtained that de- 
pends on the turbulence structure. To be more 
specific, assume that the dependence of R^ on the 

beam separation £ and the time separation r might be 
approximated through a separation of variables. 

Within the validity of this assumption, all turbulence 
parameters can be obtained which are commonly 
derived from two point product mean values [ 3] . 

This is demonstrated in graphical rather than 
mathematical form. 

The integral scale of turbulence follows from the 
space correlation as shown in Figure 2. It is equal 
to that separation distance which makes the cross- 
hatched areas equally big. The integral scale depends 
on the direction of the beam separation. The sub- 
script x denotes a separation in the streamwise di- 
rection. If the beam separation is repeated along the 
remaining perpendicular y and z axis, the values L^ 

and L could be obtained. These values are used to 
z 

calculate mean- square values from the measurements 
at zero-beam separation. 
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FIGURE 2. TURBULENCE SCALES AND 
INTENSITIES FROM CROSSED-BEAM COVARIANCE 
MEASUREMENTS 

Using the existing computer programs [ 4] for 
random vibration analysis, a time log can be intro- 
duced between the records of the two photo-detector 
outputs. PLotting the calculated G values against 
this time separation instead of space separation gives 
additional turbulence parameters as shown in Figure 
3. For zero space separation, the measured, tem- 
perature-correlation function resembles the auto 
correlation and can be used to calculate the shape 
of the power spectrum. Plotting the temporal 
correlation function G (t ) £ for one constant space 
separation enables the bulk convection speed U c to 

be read from the time lag corresponding to the first 
maximum correlation. Repeating this plot for 
several beam separations gives an envelope which 
indicates the eddy-lift time. 

The pointwise estimates are only approximately 
correct because of the assumed separation of vari- 
ables. The case, which is not covered, occurs when 
the separation of variables is erroneous for turbulent 
fields with large scales. Large lateral scales now 
indicate that only relatively few wave -number 



FIGURE 3. LOCAL POWER SPECTRA 
CONVECTION SPEEDS AND EDDY LIFETIMES 
FROM CROSSED-BEAM CORRELATION 
MEASUREMENTS 

components are involved in the turbulent exchange. 
Separation of variables assumes that the lateral- 
wave-number components do not interact. Thus, the 
case not covered would demand a strong interaction 
between a few lateral-wave-number components which 
are not communicated to other regions of the wave- 
number space. Such a process will be broken up by 
the action of the nonlinear terms in the equations of 
motion. Therefore, the crossed-beam methods 
should give a good approximation of point measurements 
in almost all turbulent fields of practical interest. 

IV. EXPERIMENTAL RESULTS 

The most revealing and instructive test of the 
crossed-beam concept is to compare the optical ap- 
proximation of point measurements with known hot- 
wire measurements. All measurements were taken 
in a subsonic (M=0. 2) air jet exhausting through a 
1-inch (2. 5 cm) diameter nozzle into the atmosphere 
( Fig. 4) . Light extinction was achieved by spraying 
a small amount of liquid nitrogen into the settling 
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chamber. This produces small water droplets in the 
exhaust flow which attenuate the crossedbeams by 
scattering. 
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Jarrel Ash 
UV 

Monochromator 




FIGURE 4. THE FIRST CROSSED-BEAM 
CORRELATOR 

A diagram of the optical and electronic hardware 
is shown in Figure 4. A plane and a spherical mirror 
collects the light from a powerful mercury arc and 
transmits it as a parallel beam through the jet. A 
similar mirror combination projects the arc image 
on the entrance slit of a grating spectrograph. The 
aperture stop in front of the spectrograph is set at 
a beam diameter of 1 mm. The photo- multiplier 
then scans the first order spectrum of the grating. 

In this way, the wavelength interval AX and the wave- 
length X are adjusted by the width of the monochro- 
mator slit and the rotation of the grating respectively. 
By exchanging light sources, gratings, and photo 
detectors, the system is able to cover the spectrum 
from the vacuum ultraviolet (X = 1200A) to the in- 
frared (X = 25jn) . 

The above light source and detector arrangement 
is used to generate two beams, one in a horizontal 
direction and one in a vertical direction (Fig. 5). 



FIGURE 5. LIGHT SOURCE AND DETECTOR 
ARRANGEMENT 


A lathe bed provides a mechanical support which 
allows the beams to be aligned parallel to the nozzle- 
exit plane. A beam separation in the streamwise 
direction gives longitudinal turbulence scales and 
convection speeds. Transport of the whole system 
allows us to repeat the measurements for all axial 
and radial traverses of the jet. However, beam 
orientation at an angle to the nozzle-exit plane is not 
possible and the associated-lateral scales have not 
yet been measured. 

The contribution of correlated, light-source 
fluctuation, created by power supply main ripple, 
was eliminated using electronic filters. Figure 6 
shows measured, temporal, cross correlation for 
various streamwise beam separations that were 
taken along the center of the shear layer 3 diameters 
downstream from the nozzle exit. Each curve is 
normalized by the maximum value of the correlation 
coefficient observed at that particular separation. 

This form of presentation was necessitated by a cer- 
tain inconsistency of peak values which are not fully 
understood. However, the peaks of the various 
curves presented clearly indicate the presence 
of a convected, turbulent pattern. Plotting the 
known space separation against the time lags that 
are indicated by the maximum temporal correlation 
function results in a line. The slope of this line in- 
dicates the convection speed as discussed previously. 


29 






F. R. KRAUSE, M. J. FISHER, AND R. E. L\RSON 



FIGURE 6. CONVECTION VELOCITY FROM 
CROSSED-BEAM CORRELATION 

Figure 7 shows the temporal, cross-correlation 
results. The only difference from the previous re- 
sults is that the beam interaction point has been moved 
radially outwards by 0.2 inch (0. 51 cm). These figures 
yield two points on the convection speed profile which is 
established from the hot-wire measurements of Davis 
and Fisher [ 5] as shown in Figure 8. The difference 
between the crossed-beam measurements and the 
hot-wire measurements fall within the scatter of the 
hot-wire points. The spatial resolution of the crossed 
beam is good considering that the two measurements 
were only 0. 2 inches (0. 51 cm) apart. 



FIGURE 7. CONVECTION VELOCITY FROM 
CROSSED-BEAM CORRELATION 



FIGURE 8. COMPARISON OF RADIAL DISTRIBUTION 
OF CONVECTION VELOCITY 


The comparison between the optically determined 
space -correlations and the hot-wire results of 
Laurence [ 6] and Bradshaw [ 7] are shown in Figure 
9. The optical results fall in between the hot-wire 



FIGURE 9. MEASUREMENT OF EDDY SCALES 

results. Once again, this good comparison is very 
encouraging. At larger space separations, the crossed- 
beam method indicates considerable negative space - 
correlation, whereas the hot-wire correlations have 
decayed to almost zero. Many of Laurence’s auto- 
correlation measurements and Bradshaw’s space cor- 
relations also show a negative loop, but its amplitude 
is generally less than that indicated by the crossed- 
beam method. 
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The discrepancies at the large space separations 
could be explained by errors in either the hot-wire 
or the crossed-beam method. The hot wire is sub- 
jected to all velocity components; and at large sepa- 
rations, the transverse-velocity fluctuations could 
add so much noise that the correlation disappears. 
The crossed-beam method could provide erroneous 
results if the water droplets follow only the large- 
scale eddies. These eddies are limited to a range of 
small wave numbers which show a large correlation 
that disappears when all wave numbers are con- 
sidered. 

Figure 10 shows an initial measurement of the 
radial intensity profile of concentration fluctuations. 



FIGURE 10. TURBULENT INTENSITY PROFILES 

For this measurement, the beams were arranged to 
intersect in a plane 6 diameters downstream of the 
exit. The vertical beam was moved to various po- 
sitions across the flow, thereby changing the radial 
location of the intersection point. For each radial 
separation, the rms amplitude was evaluated fol- 
lowing the procedure that was already described in 
Figure 2. Since lateral scales cannot be measured 
with the present support, it was assumed to be equal 
to the longitudinal scale. In spite of this simpli- 
fication, the agreement with the hot-wire results 
from Bradshaw are again very encouraging. 

The crossed-beam results in Figure 10 are 
scattered. Any one determination shows a com- 
paratively smooth variation of intensity and the peak 
value occuring close to the center of the shear layer 
as expected. However, the decrease on either side 
of the maximum varied in repeated runs. Since the 
runs were performed on days when the external re- 
lative humidity was very different, it is suspected 


that the water vapor content of the air delivered by 
the jet modified appreciably the degree of persistence 
of water droplets in the flow. Thus, it is tentatively 
concluded that the scattering of the results is a real 
effect associated with a lack of control in the intro- 
duction of the water droplet content rather than in the 
limitation of the crossed-beam method. 


V. FUTURE APPLICATIONS 


The first contribution to the Saturn program will 
be to measure turbulence shock-wave interaction in 
the plume impingement and recirculation areas of the 
S-IC and S-n rocket clusters. Model tests at MSFC’s 
impulse base flow facility indicated extremely high 
fluctuations of wall pressures and heat-transfer rates 
as shown in Figure 11. These fluctuations create 



TIME - 2 milliseconds/cm 


TANK PRESSURE *55/* of H* 



TIME > 2 Milliseconds/cm 
168 /*Hg -TANK PRESSURE 



SI STAGE 
IB FF MODEL 



SI STAGE 
IB FF MODEL 


FIGURE 11. TURBULENT FLUCTUATIONS OF 
BASE PRESSURE AND HEAT TRANSFER RATES 

problems when predicting heat-transfer rates. 
Turbulence levels of only 3 percent have sometimes 
increased the stagnation point heat transfer on cyl- 
inders up to 80 percent. Hot-wire measurements on 
the two-dimensional -base -flow model already indicate 
turbulence levels in excess of 100 percent as shown 
in Figure 12. Therefore, the turbulence effects on 
heat transfer might play a dominant role. They 
could increase or decrease the heat transfer by either 
transporting the hot particles into the wall layers or 
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throwing them back from the wall. These effects 
cannot be included in the present model tests on 
analytical solutions. Therefore, the crossed-beam 
method is needed to study the turbulent heat and mass 
transport in base flows. 


Other optical fluctuation measurements in 
clustered rocket exhausts are necessary to predict 
and/or prevent high acoustical loads on the upper 
stages. The presence of powerful Mach-wave and 
shock- wave sound radiation is shown on the shadow- 
graph in Figure 13. In the case of single jet, the 
sound radiations are not so dangerous, since they 
are directed away from the vehicle. However, on a 
rocket cluster, the sound radiation can pass through 
the subsonic hole between the plume impingment 
shocks and will then hit the vehicle directly. The 
presence of such a subsonic hole and the oscillation 
of the plume impingement shock are clearly indicated 
on the schlieren and interfero grams of our two- 
dimensional-base-flow model (Fig. 14). The as- 
sociated, aerodynamic, feedback loops will produce 
a jet noise problem on the upper stages, as soon as 
plume impingement and a choked-flow establish a 
finite base pressure. For the S-II stage, a con- 
servative estimate of sound pressures is 153 dB. 


FIGURE 12. RMS VELOCITY PROFILES IN 
TWO-DIMENSIONAL- BASE FLOW WITH AND 
WITHOUT WALL EFFECT 
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The use of laser beams, if feasible, would allow 
crossed-beam measurements over vast areas. 
Speculations could also be made about water vapor, 
carbon dioxide, and ozone concentrations in the lower 
and upper atmosphere. A study of intense -radiation 
belts is conceivable which use only detectors and no 
light sources. 

If the crossed-beam method works on both 
particle scattering and gaseous absorption, it could 
be used on two -phase flows. Fuel mixing in 


combustion chambers, fuel droplets, and ion particles 
in liquid and solid rocket exhausts, gas bubbles, and 
condensation in tanks, etc. , are some of the ap- 
plications that might be of great interest to other 
laboratories. 

The present development of technique is aimed 
at the oxygen absorption in the untraviolent because 
of the importance of air flows. The insufficient 
storage capacity of the computer memory is the only 
limitation to getting accurate cross -correlation 
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estimates. However, a new computer program will 
eliminate this difficulty. 

VI. CONCLUSIONS 


Local fluctuation measurements through a space 
and time correlation of optical signals have been suc- 
cessful in a small subsonic jet where the light ex- 
tinction was produced by scattering at water droplets. 
Comparing the hot-wire data with the crossed-beam 
data lead to the following conclusions: 

(a) Probe interference and cancellation problems 
were avoided. 

(b) One shot estimates of area integrated cor- 
relation functions can be made. 

(c) Approximations can be made on the local 
power spectra, convection speeds, turbulence 
scales, and eddy lifetimes. 


(d) The crossed-beam method can work with 

existing random vibration computer programs. 
An extension that includes light extinction 
by gaseous absorption and/or emission 
would show the following potential: 

(1) The adjustment to special theromo- 
dynamic properties, such as species 
concentrations by a proper selection of 
optical wave length and wavelength 
interval. 

(2) Two phase -flow applications by inde- 
pendent measurement of droplet scat- 
tering and gaseous absorption. 

(3) The measurement of combustion in- 
stability and plasma instability. 

(4) The location and intensity of atmospheric 
concentration fluctuations and radiation 
belts. 
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UNITS OF MEASURE 


In a prepared statement presented on August 5, 1965, to the 
U. S. House of Representatives Science and Astronautics Committee 
(chaired by George P. Miller of California), the position of the 
National Aeronautics and Space Administration on Units of Measure 
was stated by Dr, Alfred J. Eggers, Deputy Associate Administrator, 
Office of Advanced Research and Technology: 

"In January of this year NASA directed that the international 
system of units should be considered the preferred system of units, 
and should be employed by the research centers as the primary 
system in all reports and publications of a technical nature, except 
where such use would reduce the usefulness of the report to the 
primary recipients. During the conversion period the use of cus- 
tomary units in parentheses following the SI units is permissible, 
but the parenthetical usage of conventional units will be discontinued 
as soon as it is judged that the normal users of the reports would 
not be particularly inconvenienced by the exclusive use of SI units." 

The International System of Units (SI Units) has been adopted 
by the U. S. National Bureau of Standards (see NBS Technical News 
Bulletin, Vol. 48, No. 4, April 1964) . 

The International System of Units is defined in NASA SP-7012, 
"The International System of Units, Physical Constants, and 
Conversion Factors," which is available from the U. S. Government 
Printing Office, Washington, D. C. 20402. 

SI Units are used preferentially in this series of research re- 
ports in accordance with NASA policy and following the practice of 
the National Bureau of Standards. 




